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^ ! Abstract 

a: 

O 1 We study the possibility of obtaining large direct CP violation in the charmed 

, hadron decays D + — * p + p°(uj) — * p + ir + ir~, D + — > ir + p°(uj) — * Tr + ir + ir~, 

D° (j)p°(uj) -» 07r+vr-, D° -> ??p (w) -> ^vr+vr-, L>° -> r/pV) -» 
7]'tt + it~ , D° — > 7r°p°(u;) — ► 7r 7r + 7r _ , and A c — > pp (u>) — > p7r + 7r~ via p — u> 
mixing. The analysis is carried out in the factorization approach. The CP 
violation parameter depends on the effective parameter, iV c , which is relevant 
to hadronization dynamics of each decay channel and should be determined 
by experiments. It is found that for fixed N c the CP violation parame- 
ter reaches its maximum value when the invariant mass of the tt + tt~ pair 
is in the vicinity of the uj resonance. For D° — » ir°p (u>) — > 7r 7r + 7r~ and 
A c — > pp°(uj) — ► p7r + 7r _ the maximum CP violating asymmetry is bigger 
than 1% when 1.98 < iV c < 1.99 and 1.95 < N c < 2.02, respectively. We 
also estimate the decay branching ratios for D" — > and A c — > p/3 for 
these values of N c , which should be tested by future experimental data. 



PACS Numbers: 11.30.Er, 13.20.Fc, 14.20.Lq, 12.39.-x 



I. Introduction 

Although CP violation has been known in the neutral kaon system for more 
than three decades its dynamical origin still remains an open problem. Besides the 
kaon system, the study of CP violation in heavy quark systems has been a subject 
of intense interest and is important in understanding whether the Standard Model 
provides a correct description of this phenomenon through the Cabbibo-Kobayashi- 
Maskawa (CKM) matrix. Actually there have been many theoretical studies in the 
area of CP violation in beauty and charm systems and some experimental projects 
have been proposed^. 

Recent studies of direct CP violation in the B meson system have suggested 
that large CP- violating asymmetries should be observed in forthcoming experiments. 
However, in the charm sector, the CP violation is usually predicted to be small. A 
rough estimate of CP violation in charmed systems gives an asymmetry parameter 
which is typically smaller than 10~ 3 due to the suppression of the CKM matrix 
element s0. By introducing large final-state-interaction phases provided by nearby 
resonances, Buccella et al. predicted larger CP violation, namely a few times 10 -3 [|J. 
On the other hand, experimental measurements in some decay channels are consis- 
tent with zero asymmetry 0. 

Direct CP violation occurs through the interference of two amplitudes with dif- 
ferent weak and strong phases. The weak phase difference is determined by the 
CKM matrix elements and the strong phase is usually very uncertain. In Refs.fB], 0, 
the authors studied direct CP violation in hadronic B decays through the interfer- 
ence of tree and penguin diagrams, where p — uj mixing was used to obtain a large 
strong phase (as required for large CP violation). This mechanism was also applied 
to the hadronic decays of the heavy baryon, where even larger CP violation may 
be possible @. In the present paper we will investigate direct CP violation in the 
hadronic decays of charmed hadrons, involving the same mechanism, with the aim 
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of finding channels which may exhibit large CP asymmetry. 

Since we are considering direct CP violation, we have to consider hadronic ma- 
trix elements for both tree and penguin diagrams which are controlled by the effects 
of nonperturbative QCD and hence are uncertain. In our discussions we will use 
the factorization approximation so that one of the currents in the nonleptonic decay 
Hamiltonian is factorized out and generates a meson. Thus the decay amplitude of 
the two body nonleptonic decay becomes the product of two matrix elements, one 
related to the decay constant of the factorized meson and the other to the weak 
transition matrix element between two hadrons. There have been some discussions 
of the plausibility of factorization P| , and this approach may be a good approx- 
imation in energetic decays. In some recent work corrections to the factorization 
approximation have also been considered by introducing some phenomenological 
nonfactorizable parameters which depend on the specific decay channels and should 
be determined by experimental datafLl|, HH HH B - 

The effective Hamiltonian for the AS* = 1, weak, nonleptonic decays has been 



discussed in detail in Refs. ||15||||16|| , where the Wilson coefficients for the tree and 
penguin operators were obtained to the next-to-leading order QCD and QED cor- 
rections by calculating the 10 x 10, two-loop, anomalous dimension matrix. The 
dependence of the Wilson coefficients on renormalization scheme, gauge and infra- 
red cutoff was also discussed. The formalism can be extended to the charmed hadron 
nonleptonic decays in a straightforward way. 

The remainder of this paper is organized as follows. In Section II we calculate 
the six Wilson coefficients of tree and QCD penguin operators to the next-to-leading 



order QCD corrections by applying the results of Refs. ||T5[ JT(| . Then in Section III we 
give the formalism for the CP-violating asymmetry in charmed hadron nonleptonic 
decays and show numerical results. Finally, Section IV is reserved for a summary 
and some discussion. 
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II. The effective Hamiltonian for nonleptonic charmed hadron decays 



In order to calculate direct CP violation in nonleptonic, charmed hadron decays 
we use the following effective weak Hamiltonian, which is Cabibbo first-forbidden, 
based on the operator product expansion, 

G 6 



n AC= i = ^j=[ £ v uq v: q { Cl o\ + C2 oi) - v ub v; b £ c ,o 4 ] + h.c. (i) 

q= d,s i=3 



Here q (i = 1,...,6) are the Wilson coefficients and the operators Oj have the 
following expressions: 

Of = Uoa^il - 75)g/3g/37 M (l - 75)c Q , 0\ = M7 M (1 - 75)^7^(1 - 7s)c, 
3 = ^(1-75)0^^(1-75)?', 4 = M a 7/,(l-75)c /3 £^7 Ai (l-75)ga, 
5 = u 7m (1 - 75 )c£gV(l +75)9', 6 = M a 7 Ai (l-7 5 )c /3 £^7 /x (l + 75)g;, 

(2) 

where a and (3 are color indices, and q' = u, d, s. In Eq.(^) 0\ and 2 are the tree 
operators, while O3 — 06 are QCD penguin operators. In the Hamiltonian we have 
omitted the operators associated with electroweak penguin diagrams. 

The Wilson coefficients, q (i = 1, 6), are calculable in perturbation theory by 
using the renormalization group. The solution has the following form, 

C(fi) = U(fi,m w )C(m w ), (3) 

where U mw) describes the QCD evolution which sums the logarithms (a s ln(m^//x 2 )) n 
(leading-log approximation) and a s (a s ln(m^ / //i 2 )) n (next-to-leading order) . In Refs. [|15[ 
p6| it was shown that U(rrii, m<i) can be written as 

U{mum2) = J, + u° {mum2) _ ^j) , (4 ) 

where U°(mi, m<i) is the evolution matrix in the leading-log approximation and the 
matrix J summarizes the next-to-leading order corrections to this evolution. 



The evolution matrices £/ c '(mi, m?) and J can be obtained by calculating the 
appropriate one- and two- loop diagrams respectively. The initial conditions, C(mw), 
are determined by matching the full theory and the effective theory at the scale rriw 
At the scale m c the Wilson coefficients are given by 

C(m c ) = U i (m c ,m b )M(m b )U b (m b ,m w )C(m w ), (5) 

where Uf (011,012) is the evolution matrix from m 2 to mi with / active flavors and 
M(m b ) is the quark-threshold matching matrix at m b . Since the strong interaction 
is independent of quark flavors, the matrices U±(m c ,oii), U^{m b ,mw) and M(m&) 
are the same as those in b-decays. Hence, using the expressions for {7° (mi, 7712), J 



and M{m b ) given in [15|[16], we can obtain C(m c ). 

In general, the Wilson coefficients depend on the renormalization scheme. The 
scheme- independent Wilson coefficients C(/i) are introduced by the following equa- 
tion 

C(/a)=(l + ^ T )cGu), (6) 

where R is the renormalization matrix associated with the four-quark operators 
Oi(i = 1, ...,6) in Eq.(||), at the scale m w . The scheme- independent Wilson coeffi- 



cients have been used in the literature J|, [17], |18| . However, since R depends on the 



infra-red regulator [|15[], C(/z) also carries such a dependence. In the present paper 
we have chosen to use the scheme-independent Wilson coefficients. 

From Eqs.d^)(0) and the expressions for the matrices Uf (mi, 777-2), Af(m&) and 
R in ||15|| ||16|| we obtain the following scheme-independent Wilson coefficients for 
c-decays at the scale m c = 1.35GeV: 

Ci = -0.6941, c 2 = 1.3777, c 3 = 0.0652, 

c 4 = -0.0627, c 5 = 0.0206, c 6 = -0.1355. (7) 

In obtaining Eq.(^) we have taken a s (mz) = 0.118 which leads to Aqq D = 0.226GeV 
and Aqq D = 0.329GeV. To be consistent, the matrix elements of the operators Oi 



should also be renormalized to the one-loop order since we are working to the next-to- 
leading order for the Wilson coefficients. This results in effective Wilson coefficients, 
4, which satisfy the constraint 

c i (m c )(0 J (m c )) = c:(O i ) tree , (8) 

where (Oj(m c )) are the matrix elements, renormalized to the one- loop order. The 
relations between d and d read 1171, ITS 



c[ — Ci, C 2 — C2, C3 — C3 — P s /3, 

4 = c 4 + P s , c' 5 = c 5 - P s /3, Cg = c 6 + P s , (9) 



where 



with 



(a s (m c )/87r)(fO/9 + G(m,m c ,g 2 ))c 2 , 

m 2 — x(l — x)g 2 



G(m, m c , q 2 ) = 4 / dxx(l — a;)ln 
Jo 



m 2 
Here g 2 is the momentum transfer of the gluon in the penguin diagram and m is 

the mass of the quark in the loop of the penguin diagram []. G(m, m c , q 2 ) has the 

following explicit expression 

ReG - 1 ' 5 

ImG = ~tt[ l + 2^-]^/l-4^. (10) 

Based on simple arguments at the quark level, the value of q 2 is chosen in the 
range 0.3 < q 2 /m 2 < 0.5 1| 0. From Eqs.(^), @ and ( pT)|) we can obtain numerical 
values of c[. When q 2 /m 2 = 0.3, 

c[ = -0.6941, c' 2 = 1.3777, 

c' 3 = 0.07226 + 0.01472i, c' A = -0.08388 - 0.04417i, 

4 = 0.02766 + 0.01472i, 4 = -0.1567 - 0.04417z, (11) 




m could be or m s . However, the numerical values of c! i change by at most 4% when we 
change m from to m s . Therefore, we ignore this difference in our calculations, setting m — m s . 



and when q /m c = 0.5, 



0.6941, c' 2 = 1.3777, 



0.06926 + 0.01483z, c' 4 = -0.07488 - 0.04448i, 



c' 5 = 0.02466 + 0.01483i, c' 6 = -0.1477 - 0.044481 



(12) 



In calculating the matrix elements of the Hamiltonian (|1|), we can then simply 
use the effective Wilson coefficients in Eqs.flTTp (|T2"D to multiply the tree- level matrix 



elements of the operators Oi(i 



6). 



III. CP violation in charmed hadron decays 



III.l Formalism for CP violation in charmed hadron decays 

The formalism for CP violation in B and hadronic decays || [7], can be 
generalized to the case of charmed hadrons in a straightforward manner. Let H c 
denote a charmed hadron which could be D^, D°, or A c . The amplitude, A, for the 
decay H c — > /7r + 7r~ (/ is a decay product) is: 



A = (n + n-f\H T \H c ) + (n+n- f\H p \H c ), 



(13) 



where H, T and H, p are the Hamiltonians for the tree and penguin operators, respec- 
tively. 

The relative magnitude and phases of these two diagrams are defined as follows: 



A = (n+n-flH^m 
A = (n+n-f\H T \H c ) 



1 + re i5 e^ 



1 + re l0 e 



i8 -%<$> 



(14) 



where 5 and <ft are strong and weak phases, respectively. <fr arises from the CP- 
violating phase in the CKM matrix, and it is axg[V u bV* b / (V uq V*)} for the c — » q 
transition (q — d or s). The parameter r is defined as 

(n+n-f\H p \H c ) 



(7r+ir-f\H T \H c ) 



(15) 
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The CP-violating asymmetry, a, can be written as: 

— 2r sin 5 sin < 



\A\ 


2 _ 


A\ 2 


\A 


| 2 +|A| 2 



1 + 2r cos 5 cos + r 2 



(16) 



It can be seen from Eq.fllBD that both weak and strong phases are needed to produce 
CP violation. Since in r there is strong suppression from the ratio of the CKM 
matrix elements, [V^bV^/(V« g V^,*)], which is of the order 1CT 3 (for both q = d and 
q = s this suppression is 0.62 x 10~ 3 , see Eqs.(|32"D and fl43| ) in III. 2), usually the CP 
violation in charmed hadron decays is predicted to be small. 

The weak phase for a specific physical process is fixed. In order to obtain 
possible large CP violation, we need some mechanism to produce either large sin 5 
or large r. p — u mixing has the dual advantages that the strong phase difference 
is large (passing through 90° at the u resonance) and well known. In this scenario 
one has 0, § 

(n + n-f\H T \H c ) = JS-U^ + ^ (17) 
(n+n-f\H p \H c ) = JH-n^pv + &p p , (18) 

where ty (V=p or u) is the tree and py is the penguin amplitude for producing a 
vector meson, V, by H c — > fV; g p is the coupling for p° — > 7r + 7r _ ; U puJ is the effective 
p — uj mixing amplitude and Sy 1 is the propagator of V, sy = s — m v + imyTy, with 
t/s being the invariant mass of the ir + n~ pair. The numerical values for the p — uo 
mixing parameter areg [5(J |21[: Refl pw (m 2 J ) = -3500 ± 300MeV 2 , Imfi pw (m 2 ) = 
-300±300MeV 2 . The direct coupling uj — > 7r + 7r is effectively absorbed into n p ^|[21|. 
Defining 

^ = r / e ^+^) ; hL = ae *«, 2lL = p e iS0 7 (19) 

where S a , 5p and 5 q are strong phases, one has the following expression for r and 5, 

re is = r , e iS q n ^ +_ P elSf>s J . (20) 
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It will be shown that in the factorization approach, for all the decay processes 
H c —> fn + n~ we are considering, ae l5a is real (see III. 2 for details). Therefore, we 
let 

ae iS " = g, (21) 

where g is a real parameter. Letting 

p e iSp = b + dj r / e i<5 9 = d + e i, (22) 

and using Eq. fl2"0|) , we obtain the following result when y/s ~ m u , 

m C + Di 

re lS = = - = , (23) 

(s - m 2 + gRelTpJ 2 + (^Imn^ + m^r J 2 

where 

C = (s-ml + gRefl^idlReflpu + b(s - ml) - cm u Yj\ 
-epmll^ + fem^ + c(s - ml)}} 
+ (glmfl puj + m ul T ul ){e\ReU. pu + b(s - ml) - cmjj 
+G?[Imn pw + &m w r w + c(s - ml)]}, 

D = (s -ml + gRefl p u){e\Re1l p u + b(s -ml) - anj:^} 
+d[lmR pu] + bm^Y^ + c(s - ml)}} 
-(glmfl^ + m w r a; ){(i[ReII pti ; + b(s - ml) - cmjj 
-eilmll^ + bm^ + cis-ml)}}. (24) 

The weak phase comes from [T4&V^/ (V uq V* q )}. If the operators Of, 0$ contribute 
to the decay processes we have 

sm0L = 



[p + A 2 A 4 (p 2 + rj 2 )} 2 + T] 2 



cos0 d = =, (25) 

'[p + A 2 A 4 (p 2 + r? 2 )] 2 + r? 2 
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while if 0{ and Of contribute, we have 

V 



smc 



Vp 2 + V 2 



cos0| s = - P (26) 

VP + V 



where we have used the Wolfenstein parametrization [^2j for the CKM matrix ele- 
ments. In order to obtain rsin5, rcos5 and r we need to calculate j3e l&13 and r'e lSq . 
This will be done in the next subsection. 

III.2 CP violation in H c -> /tt+tT 

In the following we will calculate the CP- violating asymmetries in H c — ► /7r + 7r~. 
In the factorization approximation p°(u) is generated by one current which has 
the proper quantum numbers in the Hamiltonian in Eq.(|Ij). In the following we 
will consider the decay processes D + — > p + p°{uj) — > p + n + n^, D + — > n + p°(uj) — > 

ri'7T + 7T~ , D° — > 7r°p°(ct;) — > 7r°7r + 7r _ , and A c — > pp°(uj) — > £>7r + 7r~, individually. 

First we consider .D + — > p + p°(uj). After factorization, the contribution to t£ + 
(the superscript denotes the decay product / in if c — ► f-n + -n~) from the tree level 
operator Of is 

(pVW|£> + ) = (p°|(^)|0)(p + |(«c)| J D + ) = 7\, (27) 

where (oW) and (uc) denote the V-A currents. If we ignore isospin violating ef- 
fects, then the matrix element of O^ is the same as that of Of. After adding the 
contributions from Fierz transformation of Of and 0$ we have 

tf = (c[ + ^ 2 )(l + l/N c )T 1 , (28) 

where we have omitted the CKM matrix elements in the expression of t? + . Since in 
Eq.(^) we have neglected the color-octet contribution, which is nonfactorizable and 



difficult to calculate, N c should be treated as an effective parameter which depends 
on the hadronization dynamics of different decay channels. In the same way we find 
that t£ + = — , so that, from Eq. ([l9|) , we have 



(ae 



-1. 



(29) 



The penguin operator contributions, p? + and p£ + , can be evaluated in the same 
way with the aid of the Fierz identities. From Eq. (|T9"D we have 



{fie 



(30) 



and 



where 



(r'e^) p4 



(4 + c4)(l + ^) + c / 5 + ^-4 
(ci + c' 2 )(l + ^-) 



v ub v; b 



v ud v; d 



v ub v c 



cb 



v ud v; d 



2 \4 



A 2 X 



p2 _|_ ^2 



1 - A 2 /2 \ (1 + A 2 \ 4 p) 2 + A 4 AV ' 



Next we consider A c — > pp°(u). Defining 



(31) 



(32) 



(p» P \0(\K c ) = (p»\{dd)\Q)( P \{uc)\A c )=T 2 



(33) 



we have 



t p p = (c[ + w c' 2 )T 2 . 



(34) 



After evaluating and the penguin diagram contributions we obtain the follow- 
ing results, 

(ae i5 «) p = -1, 



((3e 



i5f3\p 



C 4 + N c C 3 



( 2 + + + + 2(4 + wA) 



(r'e t5 ") p 



(2 + ^)4 + (l + ] f-) C ' 4 + 2( C ' 5 + 



c 'i + W a c ' 2 



v ub v; b 



v ud v: d 



(35) 
(36) 

(37) 
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For the decay channel D° — > <pp (uj) — > <J)tt + tt the operators Of and 0| con- 
tribute to the decay matrix elements. If we define 



(p cf>\Ot\D ) = (<j>\(ss)\0){p°\(uc)\D }=T 3 , 



(38) 



we have 



and 



where 



(r'e*<)* 



# = (4 + 



(ae'M* = 1, 



(/3e^)^ = 1, 

C 3 + j^ C 4 + c 5 + j\^ C 6 



v ub v; b 



V V* 

v us v cs 



v ub v; b 



V V* 

y us * cs 



1 - A 2 /2 



(39) 

(40) 
(41) 
(42) 

(43) 



'7T + 7T" 



For the decay channels D° — > i]p°(uj) — > i]7t + tt~ and .D — > rj'p°{u) — ► 77 
things become a little complicated. It is known that 77 and 77' have both mm + rfd and 
ss components. The decay constants, f%( v i\ and /^n, defined as 



(0\u^ 5 u\ri(r]')) = if% W) p» (0 1 sj^ 5 s\rj (?/)) = ^V)P^ 



(44) 



are different. After straightforward derivations we have 



ae 



^fy{y') ^tV) c 3 Wc C 4 ~~ c 5 ~~ i^: c 6 



v* 

Vus V cs 



(45) 

(46) 
(47) 



£ U -fs 

In the derivations of Eqs.(f|5|, [|(|[|7]) we have made the approximation that V u bV* b /V uc [V* d 
—VubVcb/VusV c * s - It is noted that the minus signs associated with c' 5 and c' 6 in Eq.(^) 
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arise because 77(77') are pseudoscalar mesons. Since the imaginary part of c' 3 {c!A is 
the same as that of c' 5 (c' 6 ), 5 q is zero. This leads to the strong phase, 5, being zero, 
in combination with Eqs.(f45|, fHf) . 



The decay constants /jJLn and ftr^ were calculated phenomenologically in Ref . p3| , 
based on the assumption that the decay constants in the quark flavor basis follow 
the pattern of particle state mixing. It was found that 

ft = 78MeV, f* = -112MeV, ft = 63MeV, ft = 137MeV. (48) 

For the decay process D + — > 7r + p°(u;) — > 7r + 7r + 7r~, two kinds of matrix el- 
ement products are involved after factorization, i.e., {p (u)\(dd)\0) (tt + \(uc)\D + ) 
and (n + \(ud)\0)(p°(u)\(dc)\D + ). These two quantities cannot be related to each 
other by symmetry. Therefore, we have to evaluate them in some phenomeno- 
logical quark models and hence more uncertainties are involved. Similarly for 
D° -> ttW) -> we have to evaluate (p° (u;) | ( J<i) 1 0) (tt° | (uc) | Z)°) and 

(7r°\(dd) |0) (p V) I (uc)\D°) seperately. 

The matrix elements for D — > X and D — > X* {X and X* denote pseudoscalar 
and vector mesons, respectively) can be decomposed as|24||, 



(X\J,\D) = (p D + Px - ml p m ^ fc) F^ + ^^ k^k 2 ), (49) 

(X*\J„\D) = \ t^ pu e^p p D p x ,V{k 2 ) + 1 \e;(m D + mx^A^k 2 ) 

mo + mx* 1 

a m k . v ^-o> a ' k , ,. m 



-(Pd +Px*) l xM{k ) —2m x *k^A 3 (k 

1- mx* K 

e ■ k 



+t 1 -^2mx*k fl A (k 2 ), (50) 

where J p is the weak current, k = pn — px(x*) and e p is the polarization vector 
of X*. The form factors satisfy the relations Fi(0) = F (0), A 3 (0) = A (0) and 
Mk 2 ) = ^^Mk 2 ) ~ ^^A 2 {k 2 ). 

Using the decomposition in Eqs.(fi9]) (|50|), we have for D + — > 7r + p°(co>), 

tf = -V2m D \p p \ \U + -U) f p F x {<) + U + ^1) UMml)} , (51) 



N r A J " " P ' \ Z N r 
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where f p and f w are the decay constants of the p and 7T, respectively, and p p is the 
three momentum of the p. 

It can be shown that = — . After calculating the penguin operator contri- 
butions, we have 

(ae i5a ) n+ = -1, (52) 



{fie 



, (53) 



/ 2<5 \7T + 



;r 



v ud v; d 



(54) 



where x is defined as 
x = 2f p F 1 (m 2 ) + 



/^( m ;)iUK) 

"V ~ N r 



+ fpF^rrip) 



+f n A (ml) 4 + 2 



mlUAoiml 



+2 



p iV c (m c + m d ) (m u + m d ) 
mlf^A {ml) 



(55) 



N c (m c + m d ) (m u + m d ) _ 

We can consider the process D° — > 7r°p°(co>) — > 7r 7r + 7r~ in the same way. We find 



iS a \TT° 



fpF^mp - UMml) 
f p F x {w? p )+UA {ml) 



{fie 



_ 2mlf^A (ml) i i , J_ / \ 
3^ (m c +m d )(m u +m d ){ 6 N c 5' 



X 



(rV d «) 



x 



v ub v; b 



v ud v; d 



(56) 

, (57) 
(58) 



[f P Fx(ml)+UA (ml)](c 1 + ±c> 2 ) 

In Eqs.(|52"|)-(|5"8|) the form factors Fi(m 2 ) and A (m^) depend on the inner struc- 
ture of the hadrons. Under the nearest pole dominance assumption, the k 2 depen- 
dence of these form factors are 



Fi(k 2 ) 



1 — k 2 /m\ 



2 • 



A (k 2 ) 



h 



1 - k 2 /m 2 Ao ' 



(59) 
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where mi = 2.01GeV and rriA = L87GeVpi| and hi and h,A are given by the 
overlap integrals of the hadronic wave functions of D and .Y(X*)|P4"[ [25[ . 

Having obtained the expressions for ae l5a , (3e l5 P and r'e lSq , for different decay 
processes, we may substitute them into Eq . (f23"D to obtain (rsin5) and (rcos<5) for 
each channel. Then, in combination with with Eqs. (|25|) and (|26|) , the CP- violating 
asymmetries a can be obtained from Eq . (|T6|) . 

III. 3 Numerical results 

In the numerical calculations, we have several parameters: q 2 , N c , and the CKM 
matrix elements in the Wolfenstein parametrization. As mentioned in Section II, 
the value of q 2 is conventionally chosen to be in the range 0.3 < q 2 /m 2 < 0.5. For 
the CKM matrix elements, which should be determined from experiment, we use 
A = 0.221, r] = 0.34 and p = -0.12 as in Ref.g]. 

The value of the effective N c should also be determined by experiment. Since 
the hadronization information is included in N c , the value of N c may be different for 
different decay channels. Furthermore, since the color-octet contribution associated 
with each operator in the Hamiltonian (|TJ) can vary, the effective N c in the Fierz 
transformation for each operator may be different. In general, nonfactorizable effects 
can be absorbed into the effective parameters a| ff after the Fierz transformation, 

a 2i = C 2i + 7 w \ C 2i-l! a 2i-l = C 2i~l + / at \ C 2ii = 1) 2, 3), (60) 
\^c)2i \^c!2i-l 



where 



1 i + 6 (i = l,...,6), (61) 



{N c )i 3 

with being the nonfactorizable effects, which may be different for each operator. 
However, since we do not have enough information about the operator dependence 



of £j, we assume £j is universal for each operator |13| and hence for each operator we 
use the same effective N c (= (iV c )j). 
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In the numerical calculations, it is found that for a fixed iV c there is a maximum 
point, a max , for the CP violating parameter a, when the invariant mass of the tt + tt~ 
pair is in the vicinity of the u resonance. We have calculated a max in the range 
N c > for different decay channels. In the calculations we use the following two 
sets of form factors p4|: 



Set l:h 1 = 0.69, h Ao = 0.67, h v = 1.23, h Al = 0.78, h M = 0.92, 

Set 2:hx = 0.78, h Ao = 0.77, h v = 1.55, h Al = 0.98, h M = 1.27. 

The above two sets of parameters correspond to taking the average transverse mo- 
mentum of the constituents in the meson to be 400MeV or 500MeV, respectively |24]] . 
The k 2 dependence of Ai(k 2 ), A 2 (k 2 ) and V(k 2 ) are the same as in Eq. (|59|) . 

The numerical results show that for D + — > p + p°(uj) — > p + 7r + 7r~, in the whole 
range N c > 0, we have a max < 3 x 10~ 4 , which is small. For D° — > r]p°(uj) — > t]7t + tt~ 
and D° — > r/ 1 p Q {uj) — > 7] f 7T + 7T~ , from Eqs.p5], ^S, ^) it can be seen that the strong 



phase 5 is zero. Therefore, we do not have CP violation in these decays in our 
approach. However, for other processes there is a small range of iV c in which we 
may have large a max (> 1%). 

For D° -> <f>p°{u}) -> 07T+7T-, the range of A^ c for a max > 1% is 1.98 < A^ c < 1.99, 
while for A c — > pp°(uj) — > ptt + tt~ the range is 1.95 < N c < 2.02. For D + — > 
7i + p°(uj) — > 7r + 7r + 7r _ we find that for the first set of form factors when N c > 56, 
o-max > 1%, while for the second set of form factors when N c > 136, a max > 1%, 
in the range 0.3 < q 2 /m 2 < 0.5. For D° — > n°p°(uj) — > 7r°7r + 7r~ we find that when 
1.98 < N c < 1.99 we have a max > 1% in the range 0.3 < q 2 /m 2 < 0.5 for both sets 
of form factors. 

The above ranges for N c were obtained by the requirement that we have large 
CP violation in this range. However, whether iV c can be in this range should be 
determined by the experimental data for the branching ratio of each decay channel. 
Usually the decay rate for D — > fp° is determined primarily by the tree operators, 
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0\ and O2, which are related to v . In fact, the reason why we can find large CP 
violation in some range of N c , is that in this range v p becomes small enough so that 
r', and hence r, becomes large (see Eqs.([l9|, ^)). However, if is too small the 
decay rate it yields for D — > fp° may be smaller than the experimental data. In such 
a case, the range of N c in which we could have large CP violation will be excluded 
by the data. 

The decay widths for nonleptonic decays of D-meson can be calculated straight- 
forwardly in the quark model of Refs. f24], p5[ . Since we are considering the range for 
N c in which t* p is small, we have to take into account the penguin contributions p^ p 
as well when we calculate the decay widths. In the calculations of the decay width 
for D° — ► <f)p° we use = 237MeV. We find that for the first set of form factors 
the branching ratio is smaller than 2.3 x 10 -8 and for the second set the branching 
ratio is smaller than 3.6 x 10~ 8 in the range 1.98 < N c < 1.99. The dependence of 
the branching ratio on q 2 /m 2 is negligible. These branching ratios are much smaller 
than the experimental data (6±3) x 10~ 4 |26j which corresponds to 1.31 < N c < 1.53 
(1.41 < N c < 1.60) for the first (second) set of form factors. Similarly, when N c > 56 
the branching ratio for D + — > 7r + p° is smaller than 1.0 x 10~ 5 , while the experimental 
data is (1.05 ±0.31) x 10~ 3 |26| corresponding to 2.1 < N c < 2.9 (2.5 < N c < 3.4) for 
the first (second) set of form factors. Therefore, we cannot have large CP violation 
in D° — > (j)p°(uj) — > <f)7t + iT~ and D + — > 7i + p°(uj) — > 7r + 7r + 7r~. 

However, for the decay processes D° — > 7r°p° and A c — > pp° there are no ex- 
perimental data at present f26H . Therefore, there is still a possibility that iV c could 



be in the range required for large CP violation for D° — > n°p°(uj) — > 7r 7r + 7r or 
A c — > pp°(uj) — > pir + 7r~. The decay width for D° — > 7r°p° is calculated in the 
same way and we find that for the first set of of form factors the branching ratio is 
1.4(1.7) x 10 -8 , while for the second set the branching ratio is 1.8(2.1) x 10 -8 for 
N c = 1.98(1.99). This prediction is almost independent of q 2 /m 2 c . 

The branching ratio for A c — > pp° can be calculated with the same method as that 
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Table 1: Values of Br(D + — > p + p°) with the first (second) set of form factors and 
a max for £>+ -> -> p+vr+vr-, with g 2 /m 2 = 0.3(0.5) 





5r(£>+ -> p+p°) 


C^max 


0.5 


6.3(9.6)xl0- 2 


2.0(1. 8)xl0~ 4 


1.0 


2.5(3.8)xl0- 2 


1.4(1. 3)xl0" 4 


1.5 


1.6(2.5)xl0~ 2 


l.l(0.96)xl0~ 4 


2.0 


1.3(1.9)xl0~ 2 


8.7(7.7) xl0~ 5 


3.0 


0.95(1.5) xlO- 2 


6.1(5.3)xl0- 5 



in Ref . || , where we worked in the heavy quark limit m c — > oo and used the diquark 
model hadronic wavefunctions for both the heavy baryon, A c , and the proton, p. 
As in the neutron case, in the diquark model the Clebsch-Gordan coefficient of the 
u[ud\ component ([ud\ is the scalar diquark) is also l/y/2 for the proton P7|. We 
find that for N c = 1.95(2.02) the branching ratio for A c -> pp° is 7.9(7.5) x 10~ 9 
for b = l.nGeV- 1 , corresponding to (k]_)l = 400 MeV, and 6.9(7.0) x 10~ 9 for 
b = 1.18(jreV -1 , corresponding to (k\)^ = 600 MeV, where (k\) is the average 
transverse momentum of the c quark in the A c . Again the branching ratio is very 
insensitive to q 2 /m 2 c . 

In Tables 1-5 we list numerical results for a max and Br(H c — > fp°) for various 
processes, with different values of N c and q 2 /m 2 c . It should be noted that Br(H c — ► 
fp°) is almost same for q 2 /m 2 = 0.3 and q 2 /m 2 = 0.5. Table 1 shows explicitly 
that, for D + — ► p + p°(uj) — > p + 7r + 7r _ , a max is at most ~ 10~ 4 no matter what A^ c is. 
From Tables 2 and 3 we can see that in the region of N c allowed by the experimental 
data, a max is of the order 10~ 4 for D° — > (f)p (uj) — > (fiir + 7r~ and D + — > n + p°(uj) — > 
7r + 7r + 7r _ . It can also be seen explicitly from Tables 4 and 5 that there is a range 
for N c in which a max may be bigger than 1% for D° — > 7i°p°(uj) — ► 7r 7r + 7r _ and 

A c — > pp°(u) —>■ p7T + 7T~. 

In Fig.l we plot the numerical values of the CP- violating asymmetries, a, for 
D° -> 7r°p (w) -> 7r°7r+7r- with A" c = 1.99 and g 2 / m c = °- 3 > °- 5 ( for N c = 1-98 we 
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Table 2: Values of Br(D° -> 0p°) and a max for D° -> #V) -> ^tt+tt", with 
g 2 / m 2 = 0.3(0.5) 



First set of form factors 


iV c 


Br(D° -> 0p°) 


"max 




1.31 


9.0xl0" 4 


1.2(1.3)xl0" 


-4 


1.36 


7.2xl0" 4 


1.3(1.3)xl0" 


-4 


1.41 


5.7xl0" 4 


1.4(1.4)xl0" 


-4 


1.46 


4.4xl0" 4 


1.5(1. 5)xl0' 


-4 


1.53 


3.0xl0" 4 


1.6(1.6)xl0" 


-4 


Second set of form factors 


N c 


Br(D° -> 0p°) 


"max 




1.41 


8.9xl0" 4 


1.4(1.4)xl0' 


-4 


1.46 


6.9xl0" 4 


1.5(1.5)xl0" 


-4 


1.51 


5.3xl0" 4 


1.6(1.6)xl0' 


-4 


1.56 


4.0xl0" 4 


1.7(1.7)xl0" 


-4 


1.60 


3.0xl0~ 4 


1.8(1.8)xl0" 


-4 



Table 3: Values of Br(D + — > 7r + p°) and a max for D + — > 7r + p°(c<j) — > 7r + 7r + 7r , with 
g 2/ m 2 = o.3(0.5) 



First set of form factors 


N c 


Br(D + -> 7T+p u ) 


"max 




2.1 


1.4xl0" 3 


-3.0(-3.0)xl0" 


-4 


2.5 


9.8xl0" 4 


-3.9(-4.0)xl0" 


-4 


2.9 


7.3xl0" 4 


-4.9(-4.9)xl0" 


-4 


Second set of form factors 




Br(D + -> 7T+p°) 


"max 




2.5 


1.3xl0" 3 


-4.0(-4.0)xl0" 


-4 


3.0 


9.2xl0" 4 


-5.2(-5.1)xl0" 


-4 


3.4 


7.3xl0~ 4 


-6.1(-6.0)xl0" 


-4 
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Table 4: Values of Br(D° — > 7r°p°) with the first (second) set of form factors and 



^7T° 


p°(L0) -> 7 


r 7T + 7T- 


, with q 2 /m 2 r = 0.3(0 


.5) 


N c 


Br(D° - 


-> ttV : 


Omax (Set 1) 


Omax (Set 2) 


0.5 


2.1(2.8) 


xlO" 2 


1.0(0.94) xlO- 4 


1.1(0.93) xlO" 4 


1.0 


2.3(3.0) 


xlO" 3 


9.5(8.2)xl0- 5 


9.3(8. 1)xl0" 5 


1.5 


2.5(3.3) 


xlO- 4 


9.7(8. 7)xl0- 5 


9.3(8.4)xl0" 5 


1.9 


4.8(6.2) 


xlO" 5 


-7.1(-8.3)xl0- 4 


-7.4(-8.6)xl0- 4 


1.98 


1.4(1.8) 


xlO" 8 


-1.4(-1.6)xl0~ 2 


-1.5(-1.7)xl0" 2 


1.99 


1.7(2.1) 


xlO" 8 


1.4(1.6)xl0 -2 


1.4(1. 7)xl0" 2 


2.1 


7.3(9.4) 


xl0~ 6 


7.2(8. 1)xl0~ 4 


7.5(8.3)xl0" 4 


2.5 


1.0(1.3) 


xlO" 4 


2.5(2.8)xl0~ 4 


2.7(2.9)xl0" 4 


3.0 


2.8(3.6) 


xlO" 4 


2.0(2. 1)xl0" 4 


2.1(2. 1)xl0" 4 


10.0 


1.6(2.0) 


xlO- 3 


1.5(1.4)xl0" 4 


1.5(1.4)xl0" 4 



Table 5: Values of Br(A c -> pp°) with (kj_)v = 400 MeV (600MeV) and a max for 
A c -> W) -> pvr+TT-, with g 2 /m 2 = 0.3(0.5) 



N c 


Br(A c -> pp°) 


C^max 


0.5 


2.2(1. 9)xl0~ 4 


2.0(1.8)xl0 -4 


1.0 


2.4(2.1)xl0~ 5 


3.3(3.0) xlO' 4 


1.5 


2.6(2.3)xl0- 6 


7.3(6.7) xlO" 4 


1.9 


4.9(4.3) xlO" 8 


3.9(4.1)xl0- 3 


1.95 


7.9(6.9)xl0~ 9 


1.0(1. 1)xl0- 2 


2.02 


7.5(7.0) xlO" 9 


-l.l(-1.0)xl0 -2 


2.1 


7.5(6.5) xlO" 8 


-3.4(-3.1)xl0- 3 


2.5 


1.1(0.92) xlO -6 


-8.0(-7.4)xl0- 4 


3.0 


2.8(2.5)xl0- 6 


-4.3(-4.0)xl0- 4 


10.0 


1.6(1.4)xl0 -5 


-l.l(-l.l)xl0 -4 
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have similar results) as a function of the invariant mass of the 7t + tt" pair. It should 
be noted that in Fig.l we used the first set of form factors. The results for the second 
set change very little. In Fig.2 we plot the results for A c — > pp°(uj) — > pn + iT~, with 
N c = 2.02 and q 2 /m 2 c = 0.3, 0.5 (for N c = 1.95 we have similar results). In both of 
these plots we find that we can have a max > 1%. 

0.02 
0.015 
0.01 
a 0.005 


-0.005 
-0.01 

720 740 760 780 800 820 840 

V^(MeV) 

Figure 1: The CP- violating asymmetry for D° — > 7r°p°(c<j) — > 7r°7r + 7r~ (with N c = 
1.99) as a function of the invariant mass of the -k + it~ pair. The solid (dotted) line 
is for q 2 /m 2 c = 0.3 (0.5). 

o 

-0.002 
-0.004 
a -0.006 
-0.008 
-0.01 
-0.012 

720 740 760 780 800 820 840 

V^(MeV) 

Figure 2: The CP-violating asymmetry for A c — > pp°(uj) — > pn + 7T~ , with N c = 2.02. 
The solid (dotted) line is for q 2 /m 2 c = 0.3 (0.5). 
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IV. Summary and discussions 



The aim of the present work was to look for possibilities of large CP violation 
in charmed meson or baryon nonleptonic decays, H c — > fn + 7i^. Since CP violation 
in the charm sector is usually estimated to be very small (less than 10~ 3 ), it would 
be fascinating to find cases where the CP violation is large (> 1%). Following 
our previous work on CP violation in the b-quark system]?], EL we have studied 
direct CP violation in D + — > p + p°(uj) — > p + 7r + 7r~, D + — > 7i + p°(uj) — > 7r + 7r + 7r~, 
D° -> #V) -> 07r + 7r-, L>° -> ?7p (a;) -> riir+n-, D° -> r/pV) -> t/tt+ti-, 
Z)° — > 7r°p°(u;) — > 7r°7r + 7r _ , and A c — > pp (u) — > p7r + 7r~ via p — u mixing. The 
advantage of p — co> mixing is that the strong phase difference is large (passing 
through 90° at the to resonance), for some fixed N c . As a result, the CP- violating 
asymmetry, a, has a maximum, a max , when the invariant mass of the 7r + 7r~ pair is 
in the vicinity of the u resonance. It was found that a max depends strongly on the 
effective parameter N c . For D + — > p + p°(u) — > p + tt + tt~, D° — > r]p°(uj) — > i]7T + ii^ 
and _D° — > i]'p°(uj) — > T]'7r + 7r~, a max is small over the whole range, A" c > 0. However, 
for other processes we found that in order to have a max > 1%, N c should be in a 
particular range in which the amplitude v p becomes small enough so that we can 
have large CP violation. This is because when v is small r can become large, 
leading to large CP violation. However, whether or not N c can be in such a range 
is determined by the decay branching ratios for H c — > fp. 

The experimental data exclude the possibility of large CP violation in D + — > 
7i + p°(uj) — > 7T + 7r + 7r _ and .D — > (pp°(uj) — ► 07r + 7r~. However, since we do not have 
data for D° — > 7r°p° and A c — > pp° at present, it is still possible that we could 
have a max > 1% for D° — > 7r p°(cG>) — > 7r°7r + 7r~ and A c — > pp°(u;) — ► p7r + 7r~ via 
p — uj mixing in some small range of Af c . We estimated that in order to have large 
CP violation, the branching ratios for D° — > 7r°p° and A c — > pp° should be around 
10 -9 ~ 10 -8 (N c = 1.98 or 1.99 for D° -> 7r°p° and A" c = 1.95 or 2.02 for A c pp°). 
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It will be very interesting to look for such large CP-violating asymmetries in the 
experiments in order to get a deeper understanding of the mechanism for CP viola- 
tion. On the other hand, the smaller branching ratios will make the measurements 
more difficult. Furthermore, the study of CP violation in A c decays may provide 
insight into the baryon asymmetry phenomena required for baryogenesis. 

Our analysis can be extended straightforwardly to say, E' c — > Ap°(cj) — > A7r + 7r~, 
and also D s — > K* p°(uj) — > K*ir + n~ , if we assume SU(3) flavor symmetry. 

In the calculations of CP violating asymmetry parameters we need the Wilson 
coefficients for the tree and penguin operators at the decay scale m c . We calculated 
the six Wilson coefficients to the next-to-leading order by applying the formalism de- 



veloped in jn| [T6j and the relevant anomalous dimension matrix elements. Since we 
only considered strong penguin operators, and since the strong interaction is inde- 
pendent of flavor, the relevant formulas in [|15| , |T6[ | can be applied to c-decays directly. 
We worked with the renormalization-scheme-independent Wilson coefficients. Fur- 
thermore, to be consistent, we introduced the effective Wilson coefficients by taking 
into account the operator renormalization to the one-loop order. 

There are some uncertainties in our calculations. While discussing direct CP 
violation, we have to evaluate hadronic matrix elements where nonperturbative QCD 
effects are involved. We have worked in the factorization approximation, which 
has not been justified completely up to now. It has been pointed out that this 
approximation may be quite reliable in energetic weak decays P, |10|- There has 



also been some discussion on nonfactorizable contributions. In Ref.[Q the authors 
introduced two phenomeno logical parameters, e\ and eg, which are scale dependent 
to parametrize nonfactorizable effects. The scale dependence of e\ and e$ cancels 



that of the Wilson coefficients c\ and C2 and it leads to af 1 and a| ff . In Refs. [ H, |T5 
renormalization-scheme-independent coefficients are used and with the definition 
in Eq.(|60D an effective N c is introduced to describe nonfactorizable effects. On 



the other hand, Buras and Silvestrini|L4| demonstrated that, in the approach of 
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Ref.|T^], it is possible to find a renormalization scheme in which the nonfactorizable 
parameters e\ and eg vanish at any chosen decay scale. In principle, such a scheme 
can be determined by experimental data. However, the present data is not accurate 
enough. 

We can see from these investigations that more work is needed before we can 
judge the factorization approach. Since c-decays are less energetic than b-decays, 
we expect even more nonfactorizable effects. In the present work, as in Ref.JTl], |I3| , 
we introduced an effective value of N c in Eq.(|60|) and assumed that it is the same 
for each a| ff (i = 1, ...,6). The value of N c should be determined by experimental 
data and it will, in general, depend on the decay channel, since hadronization dy- 
namics can be different for each channel. Furthermore, its value depends on the 
Wilson coefficients to be used. We avoid the scheme dependence in Wilson coeffi- 
cients by using the scheme independent ones. However, such coefficients do depend 
on infra-red regulators and gauge. In principle, this dependence should be can- 
celed by the matrix elements of the operators. Furthermore, while discussing the 
processes D + — > n + p°(oj) — > 7r + 7r + 7r~ and D° — > 7T°p°(uj) — > 7T°7r + 7r~ we have to 
evaluate the matrix elements in some phenomenological quark model. All these fac- 
tors may lead to some uncertainty in our numerical results. However, as pointed 
out earlier, since the large CP violation we predict is mainly caused by small v in 
some range of the phenomenological parameter associated with the breakdown of 
factorization, we expect that our predictions should still provide useful guidance for 
future investigations. 
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